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We study the ground state of two-dimensional classical 
electron solids under the influence of modulation-doped im- 
purities by using a simulated annealing molecular dynamics 
method. By changing the setback distance as a parameter, 
we find that in the strong disorder limit the ground state 
configuration contains both isolated dislocations and disclina- 
tions, whereas in the weak disorder regime only dislocations 
are present. We show, via continuum elasticity theory, that 
the ground state of the lattice should be unstable against a 
proliferation of free disclinations above a critical dislocation 
density. Associated with this, the behavior of the threshold 
electric field as a function of the setback distance changes. 
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Since Wigner ||] proposed a solid phase of electrons, 
searches for a Wigner crystal (WC) have been pursued in 
various systems. In a low density region, a realization of 
the WC on a helium film has been established [^J. The 
WC is also expected to be formed in two-dimensional 
systems in very strong magnetic fields Recently a 
considerable amount of experimental study has focused 
on the search for this magnetically induced WC . All 
of these observations have been interpreted as indirect in- 
dications of a pinned WC. There is also some evidence 
of experimental observation of the WC at zero magnetic 
field in Si MOSFET @. However, disorder, which is 
ubiquitous, makes it difficult to interpret the experimen- 
tal observations. For example, the value of the depinning 
threshold electric field of the presumed magnetically in- 
duced WC for slightly different sample geometries differs 
up two orders of magnitude Q . To interpret experimen- 
tal observations, therefore, it is necessary to understand 
the pinning sources and estimate the pinning forces. One 
aspect of this problem is that random impurities pin a 
WC and at the same time deform it, introducing defects 
in the lattice. Only a few studies have been devoted to 
clarify the possible sources of pinning in the two di- 
mensional electron system confined in a heterojunction. 
Furthermore, the details of the relation between the de- 
fects generated by impurities and the pinning jq| are not 
well understood. 

In this paper, we study pinning and orientational or- 
der in a model of the classical two-dimensional electron 
system. Our calculations are most directly applicable to 
electrons on a helium film M , although some of our ba- 
sic conclusions should apply to heterojunction systems 
in strong magnetic fields, where the filling factor is quite 



small, so that quantum exchange effects are unimpor- 
tant pfl| . Disorder is introduced by modulation-doped 
donors randomly located on a plane separated by a set- 
back distance d from the two-dimensional electron plane. 
Since the charged electrons and ions are interacting by 
a 1/r potential, the effect of impurities is characterized 
by a dimensionless constant d/ao, where ao is the lattice 
constant of a perfect lattice. It is known that arbitrar- 
ily weak disorder destroys long-range translational order 
in two dimensions associated with a crystal lattice |pd| . 
However, we can study the orientational order of this 
system, and the associated change in behavior of the de- 
pinning threshold electric field. From numerical simula- 
tions with up to 3200 particles, we have observed defects 
— predominantly dislocations and disclinations — gener- 
ated by the impurity potential. In the weak disorder 
limit (large d/ao), the ground state configuration con- 
tains a quasi-long range orientational order ]12]] (ie. hcx- 
atic phase), and we observed no free disclinations in the 
system. However, as disorder increases, isolated discli- 
nations appear, destroying the quasi- long range orienta- 
tional order |L3]]. We will argue below, based on contin- 
uum elasticity theory, that above a threshold density of 
dislocations, it is always energetically favorable to create 
isolated disclinations. 

Our principal results are summarized in Figs. [I] and 
^. Fig. |l| shows the orientational correlation functions 
for different values of d/ao- One sees that the quasi- 
long range orientational order is destroyed in the strong 
disorder limit. Associated with this crossover, the be- 
havior of the threshold electric field is changed as shown 
in Fig. [2| The crossover takes place approximately at 
d/ao ~ 1.15 ±0.1, which is consistent with the vanishing 
of orientational order. To qualitatively understand this 
behavior, it is necessary to observe the motion of the elec- 
trons as they depin. (Details will be published elsewhere 
|]l4j| .) We find that the electrons tend to flow along direc- 
tions of the local bond orientation; i.e., to flow along lo- 
cal symmetry directions of the crystal. Since the system 
does not have long range orientational order, it is neces- 
sary for electron to pass through regions of great strain 
in the lattice, where the orientation changes. These re- 
gions of strain represent bottlenecks in the electron flow. 
As the disorder strength is tuned and orientational order 
changes from quasi-long range to short range, the number 
of bottlenecks proliferates and there is a sharp increase in 
the threshold field. We also note that the threshold elec- 
tric field is very sensitive to the setback distance, which 
might explain the very disparate values of this quantity 
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in experiments M. 

More explicitly, we study a system whose energy is 
given by 



E 
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where {r,} are the electron configurations, {Ri} are 
quenched donor configurations on a modulation-doped 
plane, d is the setback distance, and e is a dielectric con- 
stant. Numerically, we use a simulated annealing molec- 
ular dynamics method |l5| to find an electron configu- 
ration, {ri}, that corresponds to the ground state or, 
at least, a typical low energy metastable state. We im- 
pose periodic boundary conditions, and use the Ewald 
sum technique |l6[ to handle the long range interaction 
in computing energies and forces. Because our study fo- 
cuses in part on the depinning properties of this system, 
and it is possible that the depinning will be inhomoge- 
neous (i.e., that current paths form inside the crystal), 
we must account for the fact that electrons exiting the 
system will need to reenter it on the opposite side of the 
unit cell. To minimize any mismatch of current patterns 
at the boundaries, along the direction of the depinning 
field we juxtapose two square boxes whose impurity con- 
figurations are mirror-images of each other. 

We take e = 13 in the numerical simulations, and 
slowly lower the temperature from above the melting 
transition down to 20mK for a typical electron density 
n = 5.7 x 10 10 cm -2 . At this temperature we measure 
the orientational correlation functions, which are sum- 
marized in Fig. |l|. The correlation functions can be well 
fitted by exponential forms in the strong disorder limit 
(small d). We find that the correlation lenght £ is a 
slowly varying function for d/ao < 1; in the interval 
1.1 < d/ao < 1-2, the correlation length rises rapidly, 
suggesting a possible divergence (and an associated phase 
transition). However, once the correlation length exceeds 
our system size, fits to either an exponential form or a 
power law become possible, and it becomes difficult to 
precisely identify what value of d/ao would be the critical 
one in an infinite system. However, based on the sharp 
increase in t;(d), we suspect that it is in the vicinity of 
1.2. This estimate is also consistent with the observation 
of the structure factor in Fig. ||, where the six-fold sym- 
metry of the orientational order []l2| appears only for the 
samples with setback distance bigger than 1.1. 

We also measure the threshold depinning field. To 
do this, we shift the positions of the particles along a 
chosen direction by steps of O.Olao (lf| up to 1-3 lat- 
tice constants. After each shift 200 MD steps are taken 
to equilibrate the system, and the pinning force is ob- 
tained in the next 100 MD steps from the averaging the 
force on particles due to impurities. The threshold field 
is determined by the maximum value of the pinning force 
multiplied by the dielectric constant during the shifting 



process; error bounds may be obtained by comparing the 
heights of the several peaks that one sees. The results 
are shown in Fig. 0. We take two samples for each set- 
back distance. We also determine the threshold field for 
several samples by gradually increasing an external field 
and by identifying the point at which the currents start 
to flow. The results obtained through the two methods 
are in good agreement. The behavior of the threshold 
depinning force changes around d/ao = 1.15 ±0.1, which 
suggest that disclinations produce extra pinning of the 
lattice. 

We now describe how free disclinations might be fa- 
vored in strongly disordered samples using standard elas- 
ticity theory |l8[. Since the charged system requires 
charge neutrality at long wavelengths, we consider a 
model whose elastic energy is given by 
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where p and A are Lame coefficients, u(r) arc displace- 
ment vectors, ity = (l/2)(dui/dxj+duj/xi), and p(r) = 
Pa + a o~ 2 3p{f) 1S an effective in-plane impurity density. 
(Summation convention is assumed for repeated indices.) 
Since the longitudinal sound velocity goes to infinity for 
1 ft interaction, we take the limit A — > oo at the end of our 
calculation. This will guarantee that the electron density 
tracks the neutralizing background, which is the correct 
physics at long wavelengths. With defining a stress ten- 
sor 



n, 



2puij + XSij (V • u) , 
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we can divide Eq. (|^) into the contributions from smooth 
elastic displacements and displacements related to de- 
fects which involve singularities. 



E = E + E' , 

E ° = 7;J d 2 rSp[dp-V-v°] 
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where vP and u' represent the regular and the singular 
parts, respectively, etc. We note that all these quanti- 
ties are well-behaved as A — * oo. As with many crystal 
systems in the presence of disorder |8|,|l7]], we find that 
our electron crystal is unstable against the formation of 
highly separated dislocations even in the weak disorder 
limits |l4| ], and this is certainly the case in our simula- 
tions. As the disorder strength is increased, the density 
of dislocations increases. It is instructive to consider the 
energy of a single disclination in a complexion of disloca- 
tions. We write the defect energy as 



E' = E' n + E' f 
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where E' D , E' BD , and E' Dp are energy of the disclina- 
tion, the coupling energy between disclination and dis- 
locations, and energy of interaction of the disclination 
directly with the background impurities, respectively. In 
a finite size system with area A, we have 
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where bj and fj are Burger's vectors and the position 
vectors of the dislocations, and s = ±1 is the "charge" of 
the disclination, specifying whether it is a 5-fold of 7-fold 
defect The contribution of (E' D ), where (...) de- 

notes a disorder average, scales only as \J A In A, and so 
is negligible for large system sizes. We note that s may 
always be chosen such that E' BD < for any complexion 
of dislocations. The contribution E' BD may be estimated 
by computing the average over complexions of disclina- 
tions. Assuming for simplicity that the dislocations are 
completely uncorrelated, (b(r)-b(r f )) = (\b\) 2 8(r — r*), we 
find 
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Thus, the net energy of an isolated disclination for large 
A is 
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so it is energetically favorable on average to create a 
disclination if 



> 
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Interpreting (|6|) 2 as the density of dislocations in lat- 
tice units, we estimate the critical density of dislocations 
as rib w (l/47r)ag~ 2 above which disclinations are ener- 
getically favorable. This value is roughly in agreement 
with what we observe in our simulations, in which iso- 
lated disclinations appear when n^ 1 w 17-22 a\. We 
note, finally, that our computation of (E') essentially es- 
timates the mean value of the disclination energy, if we 
computed the distribution of E' for all possible disclina- 
tion locations. This mean scales to ±oo as A oc Nb — > oo, 
where Nb is the number of dislocations. One can show 
that for a simple Gaussian distribution of b(r), the vari- 
ance of (E 1 ) is also proportional to Nb, so that this model 
predicts a crossover behavior from a state with very few 
disclinations to one with many at the critical dislocation 
density. It is interesting to speculate that a more re- 
alistic distribution for b(r) could convert this crossover 
to a true phase transition. This would be quite consis- 
tent with our simulations, where isolated disclinations are 



observed for small d/ao, and there is an apparantly di- 
verging orientational correlation length just before they 
disappear. Details of this calculation will be presented 
elsewhere O]. 

In conclusion, we have observed a crossover between 
a weak disorder regime and a strong disorder regime in 
a study of a model Wigner crystal under the influence 
of the modulation-doped donor impurities as we tune 
the setback distance. In the weak disorder regime, we 
have a hexatic phase where a quasi-long range orienta- 
tional order is present, whereas in the strong disorder 
regime isolated disclinations destroy the orientational or- 
der. The crossover places d/ao = 1.15 ± 0.1. Associated 
with this crossover, the behavior of the threshold electric 
field changes. We argue by a continuum elasticity theory 
that this crossover can be understood as a proliferation 
of disclinations when the density of dislocations is bigger 
than a certain critical value. 
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FIG. 1. The orientational correlation functions for various 
setback distances. Different symbols represent data for sam- 
ples with different setback distance. From top to bottom, 
d/a =2.0, 1.7, 1.5, 1.3, 1.2, 1.1, 1.0, 0.9, and 0.8. 

FIG. 2. The depinning threshold field observed in various 
samples in the unit of Eo = e/(eag). The dotted lines are 
guides to the eye. 

FIG. 3. The magnitude of the structure factor, \S(G)\, in 
the reciprocal lattice vector space for samples with d/ao = 
(a)0.8, (b)1.0, (c)l.l, (d)1.2, (e)1.4, and (f)1.7. Only the 
points at which |S(G)| > ^|5 , (G)[ ma x are plotted. For large 
setback distances, a six-fold symmetry appears, indicating the 
presence of the quasi-long range orientational order. 
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Cha + Fertig Fig. 1 
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